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POLARBEAR Constraints on Cosmic Birefringence and Primordial Magnetic Fields 
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We constrain anisotropic cosmic birefringence using four-point correlations of even-parity A-mode 
and odd-parity B-mode polarization in the cosmic microwave background measurements made by 
the POLARization of the Background Radiation (Polarbear) experiment in its first season of 
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observations. We find that the anisotropic cosmic birefringence signal from any parity-violating 
processes is consistent with zero. The Faraday rotation from anisotropic cosmic birefringence can be 
compared with the equivalent quantity generated by primordial magnetic fields if they existed. The 
POLARBEAR nondetection translates into a 95% confidence level (C.L.) upper limit of 93 nanogauss 
(nG) on the amplitude of an equivalent primordial magnetic field inclusive of systematic uncertain¬ 
ties. This four-point correlation constraint on Faraday rotation is about 15 times tighter than the 
upper limit of 1380 nG inferred from constraining the contribution of Faraday rotation to two-point 
correlations of B-modes measured by Planck in 2015. Metric perturbations sourced by primordial 
magnetic fields would also contribute to the B-mode power spectrum. Using the Polarbear mea¬ 
surements of the B-mode power spectrum (two-point correlation), we set a 95% G.L. upper limit of 
3.9 nG on primordial magnetic fields assuming a fiat prior on the field amplitude. This limit is com¬ 
parable to what was found in the Planck 2015 two-point correlation analysis with both temperature 
and polarization. We perform a set of systematic error tests and find no evidence for contamination. 
This work marks the first time that anisotropic cosmic birefringence or primordial magnetic fields 
have been constrained from the ground at subdegree scales. 


I. INTRODUCTION 

The cosmic microwave background (CMB) has been 
an invaluable resource for testing fundamental physics. 
The CMB anisotropy has a polarized component that can 
be separated into even-parity (ill-mode) and odd-parity 
(B-mode) polarization [TH^. In the standard cosmolog¬ 
ical model, density perturbations at the last scattering 
surface produce both temperature and B-mode polar¬ 
ization anisotropies. Recent measurements of B-mode 
polarization with the Planck satellite [3] are consistent 
with the standard cosmological model. Density pertur¬ 
bations do not produce primordial B-mode polarization 
at first order. Generating primordial B-modes requires 
sources with parity-odd components, such as gravita¬ 
tional waves [HEIIS], cosmic birefringence (CB) O [7], 
primordial magnetic helds (PMFs) [5Ill2j. or cosmic de¬ 
fects m- Here we examine both CB and PMF physics. 

Inflation predicts a background of gravitational waves 
that would produce a primordial B-mode signal on de¬ 
gree angular scales. The amplitude of the inflationary B- 
modes is directly related to the energy scale of inflation, a 
quantity of fundamental importance for anchoring mod¬ 
els of the early Universe. On smaller angular scales, grav¬ 
itational leasing of the CMB by the large-scale structure 
along the line of sight converts B-modes into B-modes 
M- These secondary B-modes probe the mass distribu¬ 
tion of the Universe and can provide constraints on the 
sum of neutrino masses [T51 m]. In the last two years, 
a number of experiments started to measure directly the 
B-mode power spectrum, including Polarbear [T7], BI- 
CEP2 dll [IS], ACTPol [20|, and SPTpol [2l|. 

In addition to probing inflation and the large-scale 
matter distribution, precision measurements of the CMB 
B-modes promise competitive new tests for a variety 
of exotic physics. For example, B-modes constrain the 
abundance of cosmic strings and other cosmic defects 
dH nil , supersonic bulk flows |S3] , primordial magnetic 
fields [H4I2j , and parity-violating physics [S1II|. Here we 
focus on PMFs and parity-violating interactions, both of 
which lead to birefringence, i.e., a rotation of polariza¬ 
tion converting B-modes into B-modes. B-modes gen¬ 


erated by parity-violating processes can be compared to 
those generated by a PMF via Faraday rotation, thus the 
strength of the parity-violating interaction can be quan- 
tihed by an equivalent primordial magnetic field level. 
In addition, the stress energy in the PMF sources vector- 
and tensor-mode perturbations at the time of last scatter¬ 
ing, contributing to the B-mode power spectrum 

HTj. 

The discrete symmetry groups (charge conjugation, 
parity and time reversal) play an important role in the 
standard particle physics model [35]. If both charge 
and charge plus parity symmetries are not violated then 
the observed baryon-antibaryon asymmetry cannot ex¬ 
ist (35] [33]. Cosmological models which also violate 
these symmetry groups include pseudoscalar models of 
quintessence [30], which have the benefit of naturally 
explaining the smallness of the quintessence field mass 
and of its coupling to the fields of the standard parti¬ 
cle physics model. These models couple the pseudoscalar 
and electromagnetic fields; the resulting rotation converts 
B-modes into B-modes m- 

Magnetic fields exist in all gravitationally bound struc¬ 
tures in the Universe, from planets and stars to galaxies 
and galaxy clusters. Explaining the microgauss strength 
fields observed in galaxies is challenging without a pri¬ 
mordial magnetic seed field [SS] |33] coherent over a scale 
of a few megaparsecs [53]. Giving additional impetus to 
the PMF hypothesis is the claimed detection of magnetic 
fields in the intergalactic medium [531 ES]- Candidate 
mechanisms for the generation of a PMF include infla¬ 
tionary scenarios [571155] and phase transitions [53]. De¬ 
tecting a PMF would lead to important insights into fun¬ 
damental physics and the early Universe. A recent analy¬ 
sis based on the BICEP2 detection of B-mode power [T5] 
looked for evidence of PMFs in the BICEP2 B-mode 
power spectrum at degree angular scales [30]. Planck 
data limits the magnetic field strength smoothed over I 
Mpc to BiMpc < 4.4 nanogauss (nC) at the 95% con¬ 
fidence level [ 311132 ]- Comparable bounds are obtained 
from Lyman-a spectra [35] . The next generation of CMB 
polarization experiments promise order of magnitude im¬ 
provements with the ability to detect sub-nG PMFs [44l - 

137 ]. 
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In Sec. II, we use arcminute-scale CMB polarization 
data from the Polarbear experiment to constrain the 
anisotropic cosmic birefringence power spectrum. An up¬ 
per limit of an equivalent magnetic field is obtained to 
interpret this four-point correlation measurement in Sec. 
III. In Sec. IV, an upper limit on the amplitude of an 
actual PMF is also constrained by the two-point correla¬ 
tion measurement, i.e., the Polarbear B-mode power 
spectrum. The overall structure of these sections is that 
we discuss parity-violating physics, i.e., cosmic birefrin¬ 
gence, in Secs. II and III, and the primordial magnetic 
field in Sec. IV. 


II. COSMIC BIREFRINGENCE, FARADAY 
ROTATION, AND THE ROTATION ANGLE 
ESTIMATOR 

A. Birefringence and its effect on the CMB 

Cosmic birefringence - the difference in propagation of 
different polarization states - can rotate CMB polariza¬ 
tion and convert B-modes to B-modes. One proposed 
source of cosmic birefringence is a coupling between pho¬ 
tons and a pseudoscalar field </>. Such couplings arise nat¬ 
urally in modified theories of electromagnetism which in¬ 
clude a Chern-Simons term. The Chern-Simons term can 
appear in pseudoscalar models of quintessence (STJ |48] , 
with a Lagrangian 


r — ^ p 
~ 2M ' 


( 1 ) 


to nontrivial four-point correlations. In this paper, we 
use these four-point correlations to search for anisotropic 
rotations from CMB maps [^51451) . The cosmic birefrin¬ 
gence constraints in terms of an effective Faraday rota¬ 
tion are discussed at the end of Sec. III. 


B. Faraday rotation due to primordial magnetic 
fields 

The effect of the cosmic birefringence can be described 
by an equivalent PMF inducing the Faraday rotation. In 
this section, we exclusively describe the Faraday rotation 
given by a PMF. The constraints on the actual PMF will 
be described in Sec. IV. 

A PMF embedded in the photon-baryon plasma during 
recombination will Faraday rotate the plane of polariza¬ 
tion of CMB photons, providing another mechanism for 
cosmic rotation, now with a characteristic frequency de¬ 
pendence. The rotation angle along the line of sight n is 
given by [5115^ 

where t is the differential optical depth, v is the observed 
frequency of the radiation, B is the comoving magnetic 
field, e is the electron charge, and dl is the comoving 
length element along the photon trajectory. 

A statistically homogeneous, isotropic and Gaussian 
distributed stochastic magnetic field B(x) is charac¬ 
terized by a two-point correlation function in Fourier 
space [13 [S3] by 


where is the electromagnetic field strength tensor, 
and is its dual. The coupling is suppressed by a 
mass scale M. Such an interaction will rotate the linear 
polarization of the CMB by an angle [31] 


a = 



( 2 ) 


during propagation over an interval in conformal time ry, 
where (j) = d(j)/dri. This rotation of polarization of the 
CMB creates cosmic birefringence. 

If the spatial average of the field {(j)) is not zero, the ro¬ 
tation would produce nonvanishing parity-odd two-point 
(TB) and (BB) correlation functions [7]. Such correla¬ 
tions would imply the existence of a preferred orientation 
in the Universe and are not normally expected because 
of the presumed statistical isotropy of cosmological per¬ 
turbations. 

Regardless of the value of {(j)), fluctuations in the pseu¬ 
doscalar field will generate anisotropy in the rotation an¬ 
gle a, leading to a spatially varying cosmic birefringence. 
A statistically isotropic, random a(n) creates B-mode 
power [48] with an angular dependence determined by the 
rotation power spectrum. Inhomogeneous cosmic bire¬ 
fringence also correlates the B- and B-modes, leading 


(B,(k)B,(k')) = (2^)3<5(3)(k + k')[(<5,,-fc,fc,)B(fc)], (4) 

where S{k) is the symmetric magnetic field power spec¬ 
trum and ki is a normalized component of a wave vector 
k. The antisymmetric component describes the helicity 
of the magnetic field which does not contribute to the 
Faraday rotation spectrum and has a subdominant con¬ 
tribution to the CMB power spectra, hence we omit the 
antisymmetric contribution as it is inconsequential for 
the purpose of this paper [51H55] . The shape of S{k) 
depends on the mechanism responsible for production of 
PMF and generally is taken to be a power law up to a 
certain dissipation scale /cdiss- Namely, S{k) oc fc" for 
0 < fc < kdiss and zero for k > kdiss where fcdiss de¬ 
pends on the amplitude and the shape of the magnetic 
field’s spectrum and n is the spectral index. For nearly 
scale-invariant spectra that produce CMB anisotropy, the 
value of fcdiss is irrelevant and we take it to be infinite. 
To quantify the tangling scale of the PMF we smooth 
its comoving amplitude over a length A, obtaining B\. 
For scale-invariant fields, this quantity is independent of 
A and is Beff = y/SncB, where es is the total magnetic 
energy density. 

Faraday rotation (FR) happens concurrently with the 
generation of CMB polarization during recombination. 
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However, Ref. m demonstrated that FR can be ap¬ 
plied in a second step (i.e., first produce E modes and 
then Faraday rotate them by a PMF) without introduc¬ 
ing signihcant errors. In this approximation, the power 
spectrum of the FR angle can be written as m 


^aa 




12L + 1 


TLiik) 


L + 1 
2L + 1 




( 5 ) 


where A^(fc) = k^S{k)[3c'^v ^/(IbTr^e)]^ contains all 
the physics relevant to PMF, 7z,(fc) and (k) are trans¬ 
fer functions m which are independent of the magnetic 
field and only depend on the differential optical depth. 
The discussions in this section will be applied to the Fara¬ 
day rotation equivalent of the cosmic birefringence mea¬ 
surement in Secs. II and III. 


C. Quadratic estimator and previous constraints 
on rotation power spectrum 


spectrum [liiisn]- Note that if the rotation is uniform 
over the sky, it can be entirely determined by and 

[SS]. 

In this work, we focus on the anisotropic rotation 
rather than the uniform rotation discussed in Refs. m 
[551 155] . The rotation power spectrum C£“ is derived 
from a four-point correlation of E and B via [^5U5T| 

(asB(L)a^B(L')) = (27r)2j(L - L')(Cr + W 

-|- higher-order terms), (9) 

with being the Gaussian contribution to the four- 
point function [151 [SI]. 

Previous studies have focused on constraining the uni¬ 
form rotation as well as placing upper limits on degree- 
scale rotations [sniisnsij- Constraints on the anisotropic 
cosmic birefringence power spectrum have been derived 
from WMAP-7 data using (TBTB) four-point correla¬ 
tions m- In Ref. [55], the two-point real-space correla¬ 
tion function was used to probe the anisotropic rotation. 
Both of these analyses limit the anisotropic rotation an¬ 
gle on large scales to be less than a few degrees. 


A CMB polarization experiment measures Stokes pa¬ 
rameters Q and U at different points on the sky. 
Anisotropic cosmic birefringence adds a phase factor 
g± 2 ia(n) .|.p underlying primordial CMB polarization. 
The Stokes parameters transform as 


III. BOUNDS ON ANISOTROPIC ROTATION 
FROM POLARBEAR 

A. Data analysis 


(Q ± iU) (n) = (Q ± iU ), ( 6 ) 

where Q or U denotes the primordial Gaussian CMB po¬ 
larization map, Q and U are the observed Stokes parame¬ 
ters, and Q;(n) is the anisotropic rotation field. The CMB 
polarization defined in Eq. is rotation invariant and 
can be decomposed into electric- {E-) and magneticlike 
{B-) modes [T] as 

[E±iB]{l) = /dn[Q(n)±iC/(n)]e^2i0.g-*i n^ (7) 


where 01 is the angular separation between n and 1. 

Taylor expanding the rotated CMB polarization to 
first order in the rotation angle reveals that the off- 
diagonal elements of the two-point correlation functions 
of E- and R-modes are proportional to the rotation held, 
a(n). Quadratic estimators take advantage of this fea¬ 
ture to measure the anisotropic rotation [491 - 1511 [58] . The 
quadratic estimator for CMB polarization is 


CtEBi^) = Aeb{L) 






2Cf^ cos 2011 . 


(jEEr<BB 




( 8 ) 

where 1, F, and L are coordinates in Fourier space with 
L = 1-l-F. The angular separation between 1 and 1' is 0ii', 
is the theoretical primordial power spectrum, 
and are E- and R-mode power spectra that include 
experimental noise, and Aeb{L) is a normalization fac¬ 
tor to give an unbiased estimate of the rotation power 


The PoLARBEAR telescope is located in the Atacama 
Desert in northern Chile and observes in a band cen¬ 
tered at 148 GHz with a beam size of 3.5' full width at 
half maximum. This analysis uses data on three regions 
selected for their low dust emission, hereafter referred 
to as RA4.5, RA12, and RA23 based on their right as¬ 
censions m- The total area of the three patches is 25 
square degrees, and the patches were observed by the 
POLARBEAR experiment during June 2012 to June 2013. 
This data is referred to as the first-season Polarbear 
data. 

The time ordered data (TOD) from the detectors are 
filtered and coadded into maps as described by Ref. m- 
We first flag and remove data affected by spurious instru¬ 
mental or environmental effects. The TOD are bandpass 
filtered with the upper band edge set by a low-pass filter 
and the lower band edge set by the subtraction of a first- 
order polynomial from each constant-elevation, constant- 
velocity subscan. 

A ground template, fixed in azimuth, is also removed. 
Bright radio sources are masked before removing the 
ground template and polynomial. Each pixel consists of 
two bolometers sensitive to orthogonal polarization; data 
from these two bolometers are summed and differenced 
to derive temperature and polarization TOD from each 
pixel. The TOD are then coadded with inverse variance 
weighting into maps according to a weight estimated from 
the average power spectral density between 1 and 3 Hz 
of the filtered TOD. 
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We construct an apodization window from the 
smoothed inverse variance weight map. Pixels with an 
apodization window value below 1% of the peak value 
are set to zero, as are pixels within 3' of bright sources 
in the Australia Telescope 20 GHz Survey [SS]. Q and 
U maps are transformed to E and B maps using the 
pure-B transform m- The instrument polarization an¬ 
gle is calibrated using the patch-combined power 

spectrum [niEi], so the monopole contribution to the 
anisotropic rotation is removed. 

We reconstruct the rotation field by applying the es¬ 
timator in Eq. Q to the coadded Polarbear maps 
for 1 , 1 ' G {500,2700}. The reconstructed rotation power 
spectrum is calculated as follows: 

Cr = {HL)a*{L))-Ni°^)/n, ( 10 ) 

with both the Gaussian bias and the transfer 

function Ti calculated using simulations. The mean 
estimated rotation is subtracted from the reconstruc¬ 
tions and the realization-dependent Gaussian bias is sub¬ 
tracted for the final results [SSJ [SS] . 

We create simulated map realizations of the theoretical 
spectra calculated by GAME [70]. For the simulated ro¬ 
tation maps, we assume a scale-invariant power spectrum 
L(L-|-l)C£“/27r = 10“^ rad^ (0.33 deg^). In the rotation 
simulation, map pixels are multiplied by a phase factor 
following Eq. ([^ to obtain rotated polarization maps. 
We convolve each realization by the measured beam pro¬ 
file and a transfer function that accounts for the filtering 
on the time stream, and add noise based on the observed 
noise levels in the polarization maps. The finite area of 
the Polarbear fields results in a window function that 
couples to large-scale modes, biasing them at low L. We 
correct this bias by calculating a transfer function from 
the ratio of the averaged reconstructed rotation power 
spectrum to the known input for L < 400. We validate 
the rotation reconstruction by correlating the estimated 
rotation maps with the input maps whose rotation power 
spectra are known. All the spectra for all patches agree 
with the input rotation power spectra. 

B. Systematic errors and null tests 

Systematic errors can generate spurious signals which 
might mimic the signals we want to probe. Possi¬ 
ble sources of systematic error in the Polarbear B- 
mode power spectrum have been studied extensively by 
Ref. |17j . In this section, we extend the foreground mod¬ 
eling of Ref. [T7| to the cosmic birefringence signal and 
present a new systematic null test specific to the birefrin¬ 
gence analysis. 

Astrophysical foregrounds might affect measurements 
of the anisotropic cosmic birefringence. We test the im¬ 
pact of foregrounds in four ways. First, we generate 
Gaussian realizations of the foreground emission due to 
galactic dust and synchrotron, and radio and dusty galax¬ 
ies. The templates and amplitudes for each term are 


taken from the default foreground model presented by 
Ref. [17]. The foreground realizations are added to sim¬ 
ulated GMB maps, and the rotation power spectra are 
estimated. We find adding foregrounds does not bias the 
result, but does negligibly increase the uncertainty by 
0 . 6 %. 

A potential concern about this treatment is that the 
radio galaxies, which are effectively unresolved point 
sources, should be drawn from a Poisson instead of Gaus¬ 
sian distribution. We therefore create a set of Poisson 
realizations drawn from the empirically determined num¬ 
ber counts dN/dS cx 5“, where a = —2.15 [7T] and S is 
the source flux. We limit the distribution to fluxes with 
1 mjy <S< 25 mJy. The upper limit is set by the point 
source detection threshold and the lower limit is chosen 
to be sufficiently small that it has no effect on the sim¬ 
ulated power. The equivalent power of the unresolved 
point sources is 7 /rK^ at I = 3000. We assume a 5% po¬ 
larization fraction and random polarization angles. We 
propagate the polarization maps of the unresolved point 
sources and find the simulated contamination to be neg¬ 
ligible. 

Finally, we perform two tests to quantify the worst- 
case impact of polarized galactic dust. First, we change 
the dust polarization spectrum from the default model 
in Ref. [T7] to the higher, empirical model of Ref. [I7j . 
As before, we generate Gaussian realizations of the dust 
polarization in our observing patches. We find that the 
higher empirical dust model from Ref. m leads to a 
larger, but still negligible increase in the uncertainty of 
1.2%. Second, we test whether the non-Gaussianity of 
the polarized dust anisotropy is important. For this, we 
use the dust template from the Planck Sky Model [72] 
and extract anisotropic dust maps for all Polarbear 
patches. The polarized dust maps are multiplied by a 
factor of 2 since the Planck dust study indicates a larger 
polarized fraction [751171]. The rotation power spectrum 
of these dust maps shifts the best-fit result by O.lcr. As 
we only have one template, it is impossible to determine 
if this is a bias or scatter. Both tests argue that the polar¬ 
ized dust has only a small effect on the measured rotation 
power spectra. We conservatively weaken the upper limit 
on the anisotropic birefringence signal by 0.2cr to account 
for foregrounds. 

Rotation fields for different patches should be un¬ 
correlated. We use this fact to test for poten¬ 
tial contamination via the “swap-patch” null test. 
We define a swap-patch rotation power spectra 

^aa.nuii ^ i_(L)apatch 2(L)) [ZS] ■ In the ab¬ 

sence of contamination, the swap-patch power spec¬ 
trum should be consistent with zero. The results 
are shown in Fig. [^ The probability-to-exceed val¬ 
ues (PTEs) are 64%, 10%, 71% for the three combi¬ 
nations (Q;RA23(L)aRAi2(L))> (aRA23(L)aRA4.5(L)) and 
(«rai 2 (L)q:ra 4 5 ( 1 ^))- We find no evidence of contami¬ 
nation. 
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FIG. 1: Swap-patch rotation power spectra are shown for 
each of the three patches. The power spectra are calculated 
from the rotation fields on different patches and the legend 
indicates a specific combination. The data show no evidence 
for systematic contamination. 


FIG. 3: The blue histogram shows the distribution of the 
amplitude Acs from null signal simulations. The red vertical 
line corres pon ds to the best-fit amplitude that minimizes the 
in Eq. (Tl|. 
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FIG. 2: The anisotropic cosmic rotation power spectra from 
POLARBEAR ’s first-season data in three patches. The spec¬ 
trum of an individual patch is indicated by the green (RA23), 
blue (RA12) and orange (RA4.5) colors. The coadded (red) 
power spectrum is consistent with zero. 


C. Results 

In Fig. we show the anisotropic cosmic birefringence 
power spectrum reconstructed using the quadratic esti¬ 
mator from data in three different patches, as well as the 
coadded spectrum. The measurement is consistent with 
zero and we do not detect any anisotropic rotation signal 
from POLARBEAR data. 

A scale-invariant rotation power spectrum is particu¬ 
larly interesting because it could result from inflation¬ 
ary fluctuations of a massless pseudoscalar |48]. For a 
scale-invariant rotation field, L{L + l)C'L/27r « const. 
We define a dimensionless amplitude parameter Acb 
as a factor relating an arbitrary scale-invariant spec¬ 


trum to a reference spectrum, for which L{L + 

l)C'fV27r = 10-4 rad^ (0.33 deg^). In the WMAP anal¬ 
ysis |61j . a scale-invariant power spectrum with an am¬ 
plitude 6 X 10“^ rad^ (21 deg^) is adopted. The best- 
fit amplitude of the scale-invariant anisotropic rotation 
power spectrum corresponds to the minimum of 


x\Acb) = - AcBClf) 

bb' 

where b is the index of the rotation band power and 
is the measured spectrum in band b. The covariance ma¬ 
trix Mhf,! is calculated from simulations with no cosmic 
birefringence signal. The posterior distribution is shown 
in Fig. I 

An upper limit on the amplitude of the rotation spec¬ 
trum can be interpreted as a bound on the magnitude of 
FR and the magnetic field spectrum. A scale-invariant 
PMF results in a scale-invariant FR spectrum [57]. At 
the POLARBEAR frequency v = 148 GHz, the measured 
95% confidence limit Aqb <3.1 translates into a four- 
point correlation bound on the strength of an equivalent 
PMF: i?iMpc < 90 nG, according to the relation i?iMpc = 
(2.1 X 102 nG)(j^/30 GHz)2^L(L-b l)C'g“/27r |46l |47]. 
Including estimates for known systematic errors, this 
limit becomes i?iMpc < 93 nG. Our constraint from 
the cosmic birefringence power spectrum is roughly 15 
times lower than the recent 95% confidence level limit 
of i?iMpc < 1380 nG inferred from constraining the con¬ 
tribution of Faraday rotation to the Planck polarization 
power spectra [42]. Also, compared to the WMAP cos¬ 
mic birefringence measurement [61] , the amplitude of the 
rotation power spectrum from Polarbear is roughly 60 
times smaller. 
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FIG. 4: A representative B-mode polarization power spec¬ 
trum sourced by a scale-invariant PMF. Shown are the passive 
tensor mode (green), the compensated vector mode (orange), 
the gravitational leasing contribution (blue) and the combi¬ 
nations of the leasing and vector B modes (red) and all three 
components (magenta). The PMF contribution is based on 
BiMpc = 2.5 nG, n = —2.9, a^/apMP = 10®. The data points 
are from the Polarbear first-season B-mode power spec¬ 
trum. The third point is the 95% upper limit assuming the 
band power is positive. 


IV. CONSTRAINTS ON PRIMORDIAL 
MAGNETIC FIELDS FROM THE B-MODE 
POWER SPECTRUM 


The stress energy in the PMF sources vector- and 
tensor-mode perturbations in the metric leading to a 
frequency independent contribution to the CMB’s B- 
mode polarization |2Iij . This contribution is in addition 
to the frequency dependent FR signal discussed earlier. 
There are two potentially observable frequency indepen¬ 
dent contributions to the B-mode spectrum from a nearly 
scale-invariant PMF gnuzn]. One comes from the pas¬ 
sive, or uncompensated tensor mode, which is generated 
by the PMF before neutrino decoupling. As shown with 
the dash-dotted green line in Fig. the spectrum of this 
component is practically indistinguishable from the infla¬ 
tionary gravity wave signal. The amplitude of the ten¬ 
sor contribution is proportional to B)‘[^p^[ln(ay/apMF)]^, 
where is the scale factor at neutrino decoupling and 
apMF is the scale factor at which PMF was generated. 
The passive tensor mode is not constrained by the ex¬ 
isting Polarbear analysis, which only probes I > 500 
m- However, future measurements of at I < 100 
will probe the tensor contribution, although it will likely 
be degenerate with primordial gravitational waves. 

The PMF vector modes are more directly relevant to 
the current Polarbear data as shown by the dotted 
orange line in Fig. The B-mode power spectrum gen¬ 
erated by a scale-invariant PMF peaks around I ^ 1700, 


with the peak power given by 


i{i + i)cp^ 


27T 


1^1700 


2.5 X 10"^ 


T^lMpc 

nG 


/iK^ . (12) 


The vector-mode contribution is independent of opmf- 
Therefore, the PMF B-mode power spectrum can be 
characterized by three parameters: the PMF amplitude 
^iMpc, the epoch of PMF generation f3 = In (aj,/apMF)j 
and the PMF spectral index n, where we note that the pa¬ 
rameter P only affects the tensor mode. In what follows, 
we use the Polarbear B-mode power spectrum m to 
derive constraints on BiMpc? marginalizing over the other 
parameters. 


A. Data analysis 

Our theoretical B-mode model consists of lensing and 
the PMF vector B-modes. Polarbear data measured 
the B-mode power spectrum at 148 GHz m- We use 
the published Polarbear B-mode window functions, 
band power and band variances to construct the likeli¬ 
hood function. We assume a Gaussian likelihood for the 
Polarbear data and adopt the following priors on the 
PMF parameters: 0 < BiMpc < 10 nG, —2.9 < n < —1.5 
and 0 < P < 39. A larger prior upper limit on BiMpc is 
not necessary because constraints obtained in this anal¬ 
ysis are well below this bound. The upper prior on n 
is chosen because for high n, or “bluer” PMF spectra, 
most of the PMF energy is concentrated on small scales, 
with only negligible power on scales above 1 Mpc that 
are of relevance to our data. Thus, extending the range 
of n would make no difference for our constraints, unless 
we allow for an extremely strong PMF, which is ruled 
out. On the other hand, the spectral index has to be 
larger than —3 to avoid the divergence of the PMF power 
spectrum. We take into account the systematic contami¬ 
nation of the Polarbear B-mode power spectrum con¬ 
sidered in Ref. and investigate how the systematic 
uncertainties can potentially affect the PMF constraints. 


B. Results 

In Fig. [^we show the marginalized posterior distribu¬ 
tion function (PDF) of the PMF amplitude BiMpc- We 
take advantage of the detailed study of systematic uncer¬ 
tainties affecting the B-mode power spectrum in Ref. m 
to investigate the effects on the PMF constraints. The 
PDF without systematics is in blue and the shaded area 
indicates the shift of the PDF when all known sources 
of systematic error are included. The likelihood function 
peaks at BiMpc = 0, thus only the upper bound can be 
derived. It is determined by integrating the area and the 
vertical red line shows the 95% bound of 3.9 nG; system¬ 
atic errors have a negligible impact of 5%. We have 
examined the posterior distribution of the spectral index 










FIG. 5: Posterior distribution function of amplitude BiMpc of 
primordial magnetic field using Polarbear first-season 
measurement. The vertical line indicates the 95% confidence 
level upper limit at SiMpc < 3.9 nG. The shaded area is the 
variations introduced by both the systematic and multiplica¬ 
tive effects. 


n, and find as expected that the Polarbear data do not 
constrain the spectral index. The PMF amplitude con¬ 
straint from the first-season Polarbear S-mode power 
spectrum alone is comparable to the Planck 2015 limit 
of i^iMpc < 4.4 nG at 95% confidence level, where the 
Planck results include both temperature and polarization 
information |42| . 

We have assumed a flat prior on i?iMpc for the con¬ 
straint in Fig. following the usual convention in the 
literature. Note that, as expected with a limit, the prior 
choice has a substantial effect on the resulting poste¬ 
rior and inferred limits. We investigate a uniform prior 
in the space of the observationally constrained quantity: 
-^iMpc- The 95% C.L. upper limit for this case increases 
somewhat to 4.5 nG. Another prior, the Jeffrey’s prior 
which is uniform in log]^Q[i3iMpc/nG], is frequently used 
for parameters whose magnitude is unknown. However 
the posterior for the Jeffrey’s prior diverges in this case 
due to the lack of a reliable lower bound on HiMpc, a 
conclusion also reached in Ref. [77]. Alternatively, we 
can examine for what value of i?iMpc increases the 
by 4 (analagous to 2 a) relative to HiMpc = 0; we find 
this occurs at RiMpc = 4.4nG. Based on these tests, high 
PMF amplitudes (JJiMpc > 4.5 nG) are disfavored at the 
95% C.L. by the Polarbear R-mode power spectrum 
measurement. 


V. CONCLUSIONS 

We constrain the anisotropic rotation power spectrum 
from the first-season of Polarbear data using the four- 
point correlations of the CMB polarization. The ampli¬ 
tude of this spectrum is consistent with zero and we do 


not detect anisotropic cosmic birefringence effects. The 
amplitude of an equivalent PMF interpreted from the 
anisotropic rotation power spectrum is less than 90 nG 
(93 nG) at the 95% confidence level from the four-point 
correlation functions, without (with) systematic uncer¬ 
tainties included. 

We also use the first-season Polarbear R-mode 
power spectrum to constrain the magnetically induced 
vector-mode contribution to R-modes. We find that the 
PMF amplitude from the two-point correlation functions 
is less than 3.9 nG at the 95% confidence level, assuming 
a flat prior on the PMF amplitude. This limit increases to 
4.5 nG if we instead adopt a uniform prior on the PMF- 
sourced R-mode power. Neither the anisotropic rotation 
power spectrum nor the PMF constraints show evidence 
for significant systematic errors. 

Anisotropic cosmic birefringence directly probes the B- 
mode contribution created by the parity-violating physics 
as measured by the four-point correlations of the GMB 
polarization at different angular scales. On the other 
hand, the two-point correlation function, i.e., the B- 
mode power spectrum measures all curl-like polariza¬ 
tion patterns which could be introduced by different 
sources, such as primordial tensor perturbations, grav¬ 
itational leasing effects, PMF vector (and tensor) per¬ 
turbations, Faraday rotation, and parity-violating inter¬ 
actions. Thus B-mode power spectra can provide upper 
limits on the amplitude of PMF, and the four-point cor¬ 
relation measurement would potentially distinguish the 
rotation mechanisms of the parity-violating physics from 
PMF with the upcoming multifrequency CMB experi¬ 
ments. For example, future Planck polarization data 
could be used to measure anisotropic rotation power 
spectrum over the entire sky and possibly achieve a lower 
FR upper limit, but our results will be complementary 
since the two experiments probe the cosmic birefringence 
effects on different angular scales. 

CMB polarization data will be complemented by other 
measurements as well, such as observations of 7 -ray 
emission from blazars [78H5D] . Together, these windows 
on the fundamental physics of the early Universe will 
help characterize the nature of parity-violating physics 
and primordial magnetism. 
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